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ABSTRACT 

Aims. We use rest-frame UV spectroscopy to investigate the properties related to large-scale gas outflow, and to the dust extinction 
and star-formation rates of a sample of z ~ 2 star-forming galaxies from the Galaxy Mass Assembly ultradeep Spectroscopic Survey 
(GMASS). 

Methods. Dust extinction is estimated from the rest-frame UV continuum slope and used to obtain dust-corrected star-formation rates 
for the galaxies of the sample. A composite spectrum is created by averaging all the single spectra from our sample, and equivalent 
widths and centroids of the absorption lines associated with the interstellar medium are measured. Velocity offsets of these lines 
relative to the systemic velocity of the composite, obtained from photospheric stellar absorption lines and nebular emission lines, 
are then calculated. Finally, to investigate correlations between galaxy UV spectral characteristics and galaxy general properties, the 
sample is divided into two bins that are equally populated, according to each of the following galaxy properties: stellar mass, colour 
excess, star-formation rate, and morphology; then a composite spectrum for each group of galaxies is created, and both the velocity 
offsets and the equivalent widths of the interstellar absorption lines are measured. 

Results. For the entire sample, a mean value of the continuum slope < /J >= -1.11 ± 0.44 (r.m.s.) was derived. For each galaxy, dust 
extinction was calculated from the UV spectrum and used to correct the flux measured at 1500 A (rest-frame), and the corrected UV 
flux was then converted into a star-formation rate. We found our galaxies to have an average < S FR > = 52 + 48 Mq yr^^ (r.m.s.). 
A positive correlation between SFR and stellar mass was found, in agreement with other works, the logarithmic slope of the relation 
being 1.10 ±0.10. 

Low-ionization absorption lines, associated with the interstellar medium, and measured on the composite spectrum, were found to 
be blueshifted, with respect to the rest frame of the system, which indicates that there is outflowing gas with typical velocities of the 
order of ~ 100 km/s. Finally, investigating correlations between galaxy UV spectral characteristics and galaxy general properties, we 
find a possible correlation between the equivalent width of the interstellar absorption lines and SFR, stellar mass, and colour excess 
similar to that reported to hold at different redshifts. 

Key words. Galaxies: high-redshift - Galaxies: star formation- dust, extinction - ISM: jets and outflows - Ultraviolet: ISM 



1 . Introduction ICowie et"an {1996'), was later confirmed in many works focused 

on the M»/L (Fontana et al. 2004) and, especially, the specific 
■ ■ In the past two decades, great effort has been devoted to the study ^tar-formation rate at different redshifts ( Fontana et al. 20o| 
of the Universe in the redshift range 1< z < 3. Several works ipeulner et al.l|200lt iPerez-Gonzalez et alJbOOS. .Papovich et atl 



have shown that this is the epoch when a substantial fraction of 20O6; Santini^&; RodighieroH^laOlO). Another e"l 

galaxy mass assembly took place, and when there is a peak m ^^^^^ j^^^ ^^^^ ^^^^^ ^^^-^ -^-^^^^^ ^-^^^ jj^^ discovery of a 

the evolution of the star-formation rate density through cosmic i^tion of distant (1 < z < 2) massive eai'ly-type galax- 

time ^Lillv et al. 1996'; Madau et al. 1996; Dickinson et alj2003t -^^ ^ave the spectral p roperties of p assively evolving old 

iHoEkms&^BeacQffi 2006; Daddi et a l, 2007). Observations ap- ^^^^ ^j^^ ages of 1 - 4 Gyr (iCimatti et alJl2004 . The properties 

pear to be consistent with a differential evolution of galaxies, jj^^^^ j^^^^^ ■ ^ ^ passive-hke evolution after a major 

according to which the most massive galaxies (M > 10"Mo) ^^^^^^ ^qq ^^-i^^ short-lived (r ~ 0.1 - 0.3 Gyrfl 

formed their stars earher and more rapidly than the less mas- ^j^-burst episode occurred at 1.5 < z < 3. 
sive ones. This so-called downsizing scenario, first introduced by 
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A key role in the study of star-forming galaxies, in this criti- 
cal redshift range, is played by the optical spectroscopic surveys 
that have been made possible by the advent of 8-10 m class tele- 
scopes. The observed optical regime samples the rest-frame UV 
range for galaxies at redshift z ~ 2. The rest-frame UV spec- 
trum contains much information about the star-formation activ- 
ity of the galaxy, the hot young stellar population, and the phys- 
ical, chemical, and dynamical state of the interstellar medium. 
These pro perties have been extensively studied in th e local uni- 
verse (see 'Kinnev et al.| |l993l : iHeckman et al.lll998l) and, most 
recently, higher redshift data have been collected from both large 
surv eys (Shaplev et al. 2003; ■SavagHoetalJ|2004'; Wein er et al.l 
2009 : Steidel et al. 2010) and lensed si ngle objects (Pettini et alj 
2000 . 2002; Ouid er et all l2009l l2010h . The UV continuum of 
a star-forming galaxy (SFG) is dominated by radiation emitted 
by young massive O and B stars, therefore its UV luminosity 
could be used a s a direct e stimator of the ongoing star forma- 
tion in a galaxy (lKennicutt|[l998i) . However, it is known that the 
rest-frame UV wavelength range is strongly affected by dust ex- 
tinction. Although both the composition of dust grains and their 
distribution in galaxies are not yet completely understood, em- 
pirical extinction laws are typically used to remove the effects 
of dust extinction from the observed spectral e nergy distribu- 
tion of star-forming galaxies jCalzetti et al.ll2000l) . Furthermore, 
ICalzetti et al.l ^1994) showed that the shape of the continuum can 
be fairly well approximated by a power law cc A^, where F,i 
is the observed flux {erg i"' cm^^ A"') and /? is the continuum 
slope, and that there is a strong correlation between the slope 
and nebula r extinction m easured in the optical using the B aimer 
decrement. Meure r et al.l (11999) found a linear relation between 
the UV continuum slope, parametrized by /3, and dust extinction 
at 1600 A, expressed in magnitudes. The relation, calibrated for 
a sample of local starburst galaxies, has also been found to hold 
at higher redshifts (Reddy et al. 2006, 2010). 

Rest-frame UV spectra provide not only an estimate of the 
extinction-corrected star-formation rates, but also information 
about one of the principal mechanisms currently invoked by 
galaxy evolution models to quench the star formation in a galaxy, 
as well as enrich the inter-galactic medium (IGM) with met- 
als: large-scale gas outflows. The UV spectrum of a star-forming 
galaxy, in the wavelength range between 1 100 A and 2000 A, is 
characterized by strong absorption lines, which originate from 
absorption of the stellar radiation by the interstellar medium 
(ISM) surrounding the regions of intense star-formation. These 
lines come from different elements, predominantly carbon, sil- 
icon, and iron, at different ionization stages. Some of them 
are the result of the combination of various mechanisms, such 
as stellar winds, in addition to the interstellar component, and 
therefore their profiles may be quite complicated to analyze. A 
common property of these lines is that they have been found 
to be blueshifted at velocities of about ~ 100 km/s, with re- 
spect to galaxy rest-frame, in the spectra of galaxies in the lo- 
cal universe as well as at higher redshifts (iHeckman et al.|[T998t 
IShanlev et al.]l2003h . suggesting that the gas is flowing out from 
the galaxy. The equivalent width (EW) of these absorption lines 
has bee n fou nd to be correlated with several galaxy properties 
(IHeckman et al1ll998t IShanlev et al.ll2003l: iWeiner et al.ll2009h : 
in particular, the lines are stronger in galaxies characterized by 
higher SFRs, higher masses, and higher dust extinctions, but it is 
not yet clear on which property the outflow parameters are most 
strongly dependent. These absorption lines have been shown 
to be saturated, therefore they cannot be used to infer chemi- 
cal abundances (though, in some cases, lower limits to the col- 



umn density may be estimated). IHeckman et al.l (Il998l) points 
out that, in the case of a saturated line, the EW is only weakly de- 
pendent on the ionic column density, and more strongly depen- 
dent on the velocity dispersion of the absorbing gas. Therefore, 
the correlation between EW and dust extinction might imply that 
the velocity dispersion in the absorbing gas is larger in galaxies 
with more dust extinction. By analyzing the Na /1/15890, 5 896 A 
doublet in 32 far-IR selected galaxies, IHeckman et alj (l2000l) 
concluded that winds with larger velocity spreads are driven by 
galaxies with higher star-formation rates, which contain more 
dust. Shaplev et al. ( 2003) alternatively proposes that the corre- 
lation between the EW of interstellar absorption lines and dust 
extinction can be more directly explained by considering the 
covering fraction of dusty clouds. 

Different models of galactic-scale win ds in star-forming 
galaxies have been proposed over the years ("Chevalier & Cl egg I 
1985; Tenorio-Tagle et al. 1999; Strickland & Stevens 20^ 
Nat h & Silki I2009t' iStrickland et all I2009t iSteidel et all l2010h . 
but a general consensus remains to be found. It is generally ac- 
cepted that galactic-scale winds originate from the disruption of 
bubbles of hot gas, created by a combination of the effects of 
supernovae and stellar winds in actively star-forming regions. 
Galactic winds are also believed to have a multi-phase gas struc- 
ture, in principle providing observables at all wavelengths from 
X-ray to optical. The blue-shifted absorption features that can 
be found in the wavelength range between 1100 A and 2000 
A, originate in the warm phase of the outflow. The hot phase 
(10^ - 10^ K), emitting in the X-ray regime, is the "engine" of 
the wind and carries most of the energy. The exact role of the 
warm phase in the general structure of a galactic-scale outflow 
is still under debate. In some models, the expanding hot gas en- 
trains fragments of the cooler ambient gas, leading to the for- 
mation of clouds moving at lower velocities than the main flow, 
in which the UV blue-shifted absorption lines originate. In other 
models, the UV lines originate in ambient gas pushed by the ex- 
panding bubble in the pre-wind phase. Eventually, it is not yet 
clear if the gas, though flowing out from the galaxy, is actually 
able to reach the IGM or if, instead, it is trapped in the extended 
haloes and then accreted back onto the galaxy. To discriminate 
between different models, much help should be provided by X- 
ray observations. Evidence of outflowing gas has often been de- 
tected in star-forming galaxies, but only in the UV and opti- 
cal regimes (f or a comprehensive review on galactic winds, see 
IVeilleux et all (120051) ). 

In this paper, we analyze the rest-frame UV spectra of a 
sample of 74 star-forming galaxies from GMASS (Galaxy Mass 
Assembly ultra-deep Spectroscopic Survey), a survey expecially 
designed to search for and study both passive early-type galaxies 
and their star-forming progenitors at high redshift. We focus on 
determining the extinction-corrected star-formation rates from 
the UV luminosity, which is directly measured in the spectra. 
By using composite spectra, we also study the properties of the 
absorption spectral lines associated with the ISM, and investi- 
gate possible correlations with physical galaxy properties, such 
as stellar mass, SFR, and morphology. The paper is organized 
as follows: in Sect|2] we introduce the dataset and our selection 
criteria. In Sectl3] the principal spectral features of the spectra in 
our sample are described, with particular emphasis on the ones 
related to the search for large-scale outflows. In the following 
sections, we present our analysis of the sample, along with the 
results concerning the dust extinction inferred from the contin- 
uum slope, the SFR, and the UV absorption lines associated with 
the outflowing gas. Finally, we try to correlate the physical prop- 
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Fig. 1. Top panel: BzK plot for the objects in the GMASS to- 
tal sample (red open circles). The green dots represent the 131 
objects with a secure redshift determination from the GMASS 
spectroscopic sample. The black dots represent the 74 high- 
quality SFGs analyzed in this work. Arrows indicate upper lim- 
its for objects undetected in one or more bands. Bottom panel: 
zoom on the 74 SFG sample. 



erties of the galaxies in our sample (such as star-formation rate 
and stellar mass) with UV line properties. 

2. The sample 

2.1. The GMASS survey 

GMASS is an ESO VLT large program project based on data 
acquired using the FORS2 spectrograph. The project's main 
science driver is to use ultra-deep optical spectroscopy in or- 
der to measure the physical properties of galaxies at redshifts 
1.5 <z < 3. 



The first step in the definition of the sample was to de- 
fine a region of 6.8 x 6.8 arcmin^ (matching the field of view 
of FORS2) in the GOODS-South public image taken at 4.5/im 
with the Infrared Array Camera (IRAC) mounted on the Spitzer 
Space Telescope. The sample included all the sources with 
m4.5 < 23.0 (AB system). 

The choices of photometric band and limiting magnitude 
were imposed by several considerations. Among the IRAC 
bands, the 4.5yum one offers the best compromise in terms of 
sensitivity, PSF, and image quality. In addition, it samples the 
rest-frame near-infrared up to z ~ 3, thus allowing a selection 
that is more sensitive to stellar mass. This is also the band into 
which the rest-frame 1.6 fj.m peak of the stellar SEDs is red- 
shifted for z > 1.5. The limiting magnitude was mainly dictated 
by the observational constraints imposed by the spectrograph, in 
particular by the number of available slits with respect to the sur- 
face density of targets at z > 1 .4 detected in the field. Moreover, 
at 5 < 23.0 the selection is sensitive to stellar masses down 
to log(M/M0) ^ 9.8, 10.1, and 10.5 for z = 1.4, 2, and 3, re- 
spectively. This ensures that it is possible to investigate the evo- 
lution of the galaxy mass assembly for a wide range of masses. 
The photometric selection led to the construction of a catalogue 
of 1277 objects, the GMASS total sample. From this catalogue, 
221 objects with photometric redshift Zpimt > 1.4 were selected 
as targets for ultra-deep spectroscopy. The total allocated obser- 
vational time (145 hours) was distributed over six masks: three 
observed through the 300V grism (the blue masks) and three ob- 
served through the 3001 grism (the red masks). 

Among the 221 objects, 170 could actually be put in the 
masks. To fill the masks, we also observed objects not included 
in the spectroscopic target list, obtaining a final GMASS spectro- 
scopic sample of 250 unique objects. To ensure the feasibility of 
our spectroscopy, it was necessary to set magnitude limits B < 
26.0 and I < 26.0 for the blue and red masks, respectively. The 
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Fig. 2. Redshift distribution. The black shaded histogram repre- 
sents secure spectroscopic redshifts of the observed full GMASS 
sample; the red shaded histogram represents redshifts for our 
sub-sample of 74 high-quality star-forming galaxies. In the inset, 
the distribution of photometric redshifts for the whole GMASS 
total sample is also shown (empty histogram). 
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Fig. 3. Distributions, in the redshift range 1 .5 < z < 2.8, of stellar age, stellar mass and star-formation rate computed by SED fitting 
(using ^arastonl (l2005h models with a Kroupa IMF) for: the final sample of 74 objects analyzed in this paper (shaded red histogram); 
the GMASS total sample (open histogram); and objects with a secure redshift determination from the GMASS spectroscopic sample 
(shaded black histogram). Medians of the photometric catalogue and the 74 SFG sample are also shown (vertical black and red lines, 
respectively). 



assignment of each object to either category was made according 
to colors and photometric redshifts: red objects at intermediate 
redshift {{z - K) > 2.3, Zphoi ^ 2.5) were listed as primary tar- 
gets to be included in the red masks, while objects that are blue 
((z- K) < 2.3) or have UV absorption lines that have been red- 
shifted into the optical domain (Zpimi > 2.5) were allocated as 
primary targets for the blue masks. The actual wavelength cov- 
erage for a certain slit depends, apart from the grism and order 
separation filter, on its position in the mask in the dispersion 
direction. The maximum wavelength range covered for central 
sHts was 3300-6500 A for the blue masks and 6000-11000 A for 
the red masks. To be able to include enough background for sub- 
traction of the sky, minimum slit lengths were constrained to be 
9" and 8" for the red and blue masks, respectively, while the slit 



width was always 1" for both grisms. The spectral resolution of 
A/AA X! 600 was chosen because of the large wavelength cover- 
age that it provides and is, however, high enough to resolve and 
identify spectral features for redshift determination. In total, it 
was possible to determine secure spectroscopic redshifts for 131 
objects included in the spectroscopic target list, and for another 
16 galaxies we obtained at least a redshift estimate. 

The physical properties of the galaxies belonging to the en- 
tire GMASS total sample were estimated by fitting different 
sets of evolutionary population synthesis models to broad-band 
SEDs (extending from U band to the IRAC 8/im band) derived 
from photometric data. In particul ar, we use d the results ob- 
tained with the synthetic spec tra of iMarastonI (j2005), adopting 
a Kroupa IMF (lKroupall200lh and fixed solar metallicity. Star- 
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Fig. 4. Percentage of contributing objects per wavelength unit 
to the composite spectrum of the total sample (74 SFGs) (see 
Figure |5]i. Red dots indicate the wavelength of the most promi- 
nent spectral lines detectable in the spectrum. 



formation histories were parametrized by exponentials (e~) with 
e-folding timescales t of between 100 Myr and 30 Gyr, plus 
the case of constan t SFR. For the dust extinction, a Calzetti law 
dCalzetti et alj|200 0) was assumed. By means of a fitting proce- 
dure based on;^^ minimization, a best-fit model for each galaxy 
was derived, providing estimates for age, stellar mass, ext inction 
Av, and SFR dBolzonella et al.ll2000t IPozzetti et alj|2007h . 

2.2. Selection criteria and redshift refinement 

In line with the general aims of this work, we wished to select 
a representative sub-sample of SFG spectra from the GMASS 
sample, with confirmed spectroscopic redshift and high-quality 
spectra for rest-frame wavelengt hs between 1 100 A an d 2900 A, 
which is both the range used by ICalzetti et al.l (11994 *) to define 
the continuum slope and the interval where we can find the 
blueshifted absorption spectral lines that are tracers of the warm 
phase of galactic-scale outflows. Therefore, we selected all the 
objects from the GMASS spectroscopic sample that had been ob- 
served through a blue grism mask (B < 26.0), and have a cata- 
logue redshift 1.5 < Zspec < 2.8. These redshift and photomet- 
ric criteria ensured the selection of star-forming galaxies. As an 
additional check, we plotted t he selected objects on the colour- 
colour diagram developed bv lDaddi etaP (l2004 . The GMASS 
SFG sample analyzed in this work have been verified to lie in the 
star-forming region of the BzK plot, as shown in Fig[T] Since the 
spectral features related to star-formation activity may be "con- 
taminated" by the presence of nuclear activity, objects with a 
spectrum showing AGN features were excluded. In particular, 
we found one object displaying broad emission lines, which was 
classified as a type 1 AGN, and another object with narrow C 
IV /11549 A emission typical of type 2 AGN. After these pre- 
liminary selection steps, we ended up with an ensemble of 94 
SFGs. 

The next selection was based on the quality of the redshifts. 
Since the redshifts have been originally obtained exploiting dif- 



ferent measurement tecniques, we decided to homogenize them 
adopting a recursive procedure based on the cross-correlation of 
the single spectra with a template created from the sample itself 
(iTonrv & Davislll979h . Briefly, the peak of the cross-correlation 
function between the spectrum of a galaxy and a given template 
provides the redshift and velocity dispersion of the galaxy with 
respect to the template. To quantify the significance of the re- 
sults, a cross-correlation goodness parameter R is defined^ The 
composite spectrum created by averaging the spectra from the 
preliminary sample of 94 SFGs, using catalogue redshifts, was 
taken as a template. From the cross-correlation of the template 
with each single galaxy spectrum, a new set of redshifts was ob- 
tained. On average, the difference between these new redshifts 
and the catalogue ones is Az ~ 2 x 10"^. Then, the goodness pa- 
rameter R was checked: too low a value of R means that only a 
low-quality redshift determination is possible because the spec- 
trum is of low signal-to-noise ratio (S/N). After some tests, we 
decided to use R > 3 as the constraint determining that a redshift 
is reliable, and 20 objects, which did not satisfy this constraint, 
were excluded from the final high-quality sample. 

The 74 remaining spectra were used to produce a new 
template, which was iteratively cross-correlated with each sin- 
gle spectrum until a round was reached when redshifts did 
not change significantly compared to the previous round. We 
reached convergence after five rounds and the final redshifts 
were the ones subsequently used in this work. Figure|2]shows the 
redshift distribution of the final sample of 74 high-quality SFGs, 
with reliable redshifts measured as explained above, against the 
origi nal spectroscopic G MASS sample of 131 extragalactic ob- 
jects. IKuA et al.l (|2009|) found that the peak at z ~ 1.6 represents 
a galaxy overdensity, possibly a cluster under assembly, where a 
significant portion of galaxies have a passive spectrum and hence 
are not included in our SFG sample. 

The distributions of age, stellar mass, and star-formation rate 
for the final sample of 74 SFGs, with respect to the GMASS total 
sample, are also shown in Figure [3] in the redshift range 1 .5 < 
z < 2.8. We note that the adopted selection criteria, especially 
the cut in optical magnitude, led our sample to be biased towards 
younger ages and higher star-formation rates, with respect to the 
parent galaxy population. 

3. Rest-frame UV spectroscopic features 

The composite spectrum of the sample of 74 SFGs is shown in 
Figure |2l It was generated by stacking the individual spectra, 
shifted into their rest-frame using the redshifts defined in the 
previous section, and then normalized in the wavelength range 
1600-1800 A and rebinned to a dispersion of 1 A per pixefl 
Figure|4]shows the percentage of contributing objects per wave- 
length unit to the composite spectrum, with an indication of the 
most prominent spectral lines used in our analysis. 

The c ontinuum slope of th e composite spectrum, measured 
following lCalzetti et al] (11994 . isyS = -1.15 + 0.01, which is in 
good agreement with the mean of the values measured for the 74 
individual spectra in the sample: < /5 >= -1.11 + 0.44 (r.m.s.) 
(see the following data analysis section). 

- R is the ratio of the correlation peak height to the amplitude of the 
antisymmetric noise. 

^ The composite spectrum is publicly available at 
http://www.arcetri.astro.it/~cimatti/gmass/publicdata.html 

"* The dispersion of single spectra is 2.5 A in the observed frame. The 
rest-frame dispersion of the composite was set to be 1 A ~ 2.5 A / 1 +z , 
to match the observed pixel resolution. 
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Several types of emission and absorption lines are recogniz- 
able. The most prominent ones are the lines due to absorption, by 
the interstellar medium, of the radiation coming from stars, high- 
lighted with red (low-ionization) and green (high-ionization) 
colours in Figure|5] 

Information about the most prominent spectral lines detected 
in the spectrum is summarized in Table [1] 

3. 1 . Photospheric absorption stellar lines 

Photospheric stellar lines generate by the absorption of radiation 
produced in the stellar interiors by the most external gas layers. 
Being very weak, these lines are almost undetectable in the sin- 
gle spectra of our sample, because of the low S/N in the contin- 
uum (typical <5), but they can be both identified and measured 
in composite spectra, since the averaging procedure improves 
the S/N (>15) (see Figure [T4li. As discussed later in more detail, 
stellar photospheric lines can be used to define the velocity zero 
point of the system, with respect to which the velocities of the 
outflowing interstellar gas are calculated. 

In our study here, we were able to identify five stellar lines: 
C III ^1 176 A, O IV ^1343 A, S III ^1417 A, S V ^1500 A, and 
C III ^2297 A. 

3.2. Nebular emission lines 

Two semi-forbidden carbon lines are detectable in our spectra: 
C III] .11909 A and C II] .i2326 A, unresolved at our spectral 
resolution. These fines originate in nebular regions photoionized 
by radiation from massive O and B stars, hence can in principle 
be used, as well as photospheric stellar lines, to determine the 
system rest-frame velocity. 

The C III] /1 1909 A is known to be a doublet 
(^^1906.08, 1908.73 A). The flux intensity ratio of the two fines 
depends on the physical conditions of the em itting gas, in partic- 
ular the electron density (iKeenan et al.lll99 2). For the electronic 
densities of local star-forming regions the ratio varies between 
0.5 and 1.5, with 1.5 being the most typical value and therefore 
the one we imposed in the line measurements. 

3.3. Low- and high-ionization absorption lines associated 
with the ISM 

The strongest lines detectable in the UV rest-frame spectrum 
of starburst galaxies are low- and high-ionization absorption in- 
terstellar lines. The first set of lines are associated witfi neutral 
gas regions and caused by the absorption of stellar radiation by 
the interstellar medium, which ensures that they are the optimal 
tools for detecting outflows in the host galaxy. 

In particular, the strongest lines in our spectra are Si II 
A1260 A, C II ^1334 A, Si II /11526 A, Fe II ^1608 A and Al 
III ^1670 A, Fe II ^2344 A, Fe II ^2374 A, Fe II ^2382 A, 
Fe II .12586 A, Fe II .12600 A, plus tfie blend (at our spectral 
resolution) of O and Si II at .11303 A, and tfie doublet Mg II 
A/12796, 2803 A. Some otfier weaker lines, wfiicfi we identified 
as interstellar absorption lines, can be found in Table [1] 

Higfi-ionization lines predominantly trace gas witfi T > 
IQ'^K. Tfiree doublets are recognizible in our spectra, i.e.. 
Si IV .1.11393,1402 A, C IV .1.11548,1550 A, and Al III 
/1/11854, 1862 A. Figure |6] sfiows zoomed-in spectral regions of 
tfie composite spectrum, containing tfie CIV and SilV doublets. 
Tfie spectrum fias been sfiifted into its stellar rest-frame and tfie 
expected position of tfie spectral lines is indicated by tfie dotted 



vertical lines. Unlike low-ionization lines, more tfian one process 
causes tfie formation of tfie fiigfi-ionization doublets, including 
absorption of stellar radiation by tfie interstellar medium, as well 
as tfie contributions of stellar winds and pfiotospfieric stellar ab- 
sorption, but tfiese different components are difficult to distin- 
guisfi. A study of UV O-type stellar spectra bv lWalborn & Pane3 
(11984 ) qualitatively sfiowed tfiat tfie Si IV doublet profile is dom- 
inated by tfie interstellar component, because only blue giant and 
supergiant stars contribute to tfie wind component, wfiile tfie C 
IV doublet is dominated by tfie stellar wind feature, wfiose con- 
tribution to tfie line profile comes from stars belonging to all 
luminosity classes, from main-sequence to supergiant stars. Tfie 
stellar wind component sfiould display a P-Cygni profile, wfiicfi 
is cfiaracterized by redsfiifted emission and bluesfiifted absorp- 
tion. Since spectral resolution is too low to allow a decomposi- 
tion of tfie line profile in its componentes, fiigfi-ionization lines 
are poorly suited to tfie searcfi and study of large-scale outflows. 

3.4. Other notable spectral lines 

3.4.1. Lya 

Tfie Lya .11216 A spectral line originates from recombinations 
in H II regions and can be used to study tfie pfiysical conditions 
witfiin tfiem. Wfien detected in emission, it is tfie most promi- 
nent spectral feature in tfie spectra of fiigfi-redsfiift galaxies. 
Wfien seen in absorption, tfie line fias a very broad profile, of- 
ten bluesfiifted. In our SFG sample, 32 spectra extend blueward 
to include tfie Lya region: 8 of tfiem sfiow tfie line in emission, 
wfiile 24 display absorption. We measured tfie velocity offset be- 
tween tfie Lya in emission and tfie low-ionization interstellar ab- 
sortion lines in tfie six spectra wfiere botfi sets are detectable. A 
mean velocity offset AvVm-aAs ~ 566 + 27 km/s was found, con- 
sistent witfi tfie ~ 650 km/s reported by Sfiapley et al. (2003) for 
tfieir sample of z ~ 3 Lyman break galaxies. 

Two composite spectra were tfien created, one witfi tfie spec- 
tra witfi Lya in emission (see Figure |7]l and tfie otfier witfi tfie 
spectra witfi Lya in absorption, and tfieir continuum slope was 
measured following tfie procedure described in a following sec- 
tion. We found tfie emission Lya spectrum to fiave a bluer slope 
witfi respect to tfie otfier case (respectively, jSemLvw - -1.40±0.06 
and /3absLva = -1.04+0.03), tfius confirming tfie relation between 
UV continuum slope, and tfierefore dust extinction, and tfie Lya 
profile in SFGs, as afi-eady observed in LBGs at z ~ 3. 

3.4.2. He II .11640 A 

He II /11640 A is a typical feature produced by winds from 
Wolf-Rayet (WR) stars or AGN. WR stars are thougfit to rep- 
resent a stage in the evolution of O-type stars with masses 
M > 2QMq. They are characterized by strong winds with veloci- 
ties reaching thousands km/s, causing considerably mass losses. 
The strength of the broad He II .11640 A is likely to depend 
on the ratio of WR with respect to O-type stars, which is in 
turn strongly determined by the metallicity ([Schaerer & Vaccal 
119981: lMevneilI995l: iBrinchmann et al.ll2008bl) 7 at lower metal- 
licities O-type stars are less likely to evolve into WR, and tfie 
He II becomes intrinsecally weaker. From tfie EW of He II 
.11640 A measured on tfie composite spectrum of our entire sam- 
ple, EWneii = 1 .02 + 0.10 A , we derived a WR/(WR+0) ratio 
of 0.05 ± 0.008 dScfiaerer & Vacca 1998). 

Quite a good agre ement is found witfi tfie models of 
iBrincfimann et alJ (l2008bi) . For tfie case of nearly constant SFR, 
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Fig. 5. Top panel: average spectrum of the total sample (74 SFGs). Bottom panels: zoom of the composite spectrum, above and 
below 2000 A. Spectral lines of interest are labeled. Blue: absorption stellar photospheric lines; red: interstellar absorption low- 
ionization lines; green: interstellar absorption high-ionization lines; magenta: emission nebular hnes; black: emission and absorption 
lines associated with stellar winds; cyan: interstellar fine-structure emission lines. „ 
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Table 1. Spectral lines in a SFG rest-frame UV spectrum. 



Ion 


Vacuum 
wavelength (A) 


Rest-frame 
EW° (A) 


Origin 


cm 


1175.71 


1.70 ±0.12 


stellar (photospheric) 


N V 


1238.82, 1242.80 


-1.11 ±0.20 


ISM 


Sill 


1260.42 


1.47 ± 0.09 


ISM 


Sill* 


1264.74 


-0.36 ± 0.25 


Fine-structure emission (ISM) 


O I - Si 11* 


1302.17, 1304.37 


2.43 ± 0.26'' 


ISM 


Sill* 


1309.28 


-1.26 ±0.25 


Fine-structure emission (ISM) 


CII 


1334.53 


2.09 ±0.15 


ISM 


OIV 


1343.35 


0.27 ± 0.07 


stellar (photospheric) 


SilV^ 


1393.76 


2.02 ± 0.09 


ISM -1- stellar wind (P-Cygni profile) 


SilV^ 


1402.77 


1.42 ± 0.05 


ISM -1- stellar wind (P-Cygni profile) 


Si III 


1417.24 


0.72 ± 0.29 


stellar (photospheric) 


S V 


1501.76 


0.54 ± 0.13 


stellar (photospheric) 


Sill 


1526.71 


1.94 ±0.08 


ISM 


C IV*-" 


1548.20, 1550.78 


3.24 ± 0.16'' 


ISM + stellar wind (P-Cygni profile) 


Fe II 


1608.45 


1.31 ± 0.11 


ISM 


Hell 


1640.42 


-1.02 ±0.10 


stellar wind 


AlII 


1670.79 


1.54 ± 0.08 


ISM 


Ni II 


1709.60 


0.70 ± 0.35 


ISM 


NIV 


1718.55 


1.14 ±0.18 


stellar wind (P-Cygni profile) 


Ni II 


1741.55 


0.25 ± 0.41 


ISM 


Ni II 


1751.91 


0.39 ± 0.41 


ISM 


Sill 


1808.01 


0.76 ± 0.09 


ISM 


Si I 


1845.52 


0.53 ± 0.21 


ISM 


Al IIF 


1854.72 


1.20 ±0.21 


ISM + stellar wind (P-Cygni profile) 


Al IIF 


1862.79 


0.85 ±0.11 


ISM -1- stellar wind (P-Cygni profile) 


C III]* '' 


1906.68, 1908.68 


-1.02 ±0.15'' 


nebular 


Fe II 


2249.88 


0.60 ± 0.22 


ISM 


Fe II 


2260.78 


0.50 ±0.31 


ISM 


cm 


2297.58 


0.28 ± 0.24 


stellar (photospheric) 


CII] 


2326.00 


-1.24 ±0.20 


nebular 


Fe II 


2344.21 


2.45 ±0.13 


ISM 


Fe II* 


2365.66 


-0.74 ± 0.38 


Fine-structure emission (ISM) 


FelF 


2374.46 


1.87 ± 0.13 


ISM 


Feir 


2382.76 


2.22 ± 0.15 


ISM 


pp» TT* 
re il 




1 1 -u n 1 s 

-1 . IZ ± U. lo 


Fine-structure emission (ISM) 


Fe IF 


2586.65 


2.65 ± 0.17 


ISM'' 


Fe IF 


2600.17 


2.77 ± 0.14 


ISM'' 


Fe II* 


2612.65 


-1.07 ± 1.35 


Fine-structure emission (ISM) 


Fell* 


2626.45 


-1.17 ±0.41 


Fine-structure emission (ISM) 


MglF 


2796.35 


2.43 ±0.10 


ISM^ 


MglF 


2803.53 


2.07 ± 0.10 


ISM' 


Mgl 


2852.96 


1.75 ± 0.68 


ISM 



Notes. E'Ws measured in the composite spectrum of 74 SFGs. Tha values are negative for emission lines and positive for absorption lines. 
**' Blended at our spectral resolution. 
<''> Doublet. 

*''* The measure of the EW refers to the blend of the two lines. 
Possible stellar wind contribution. 
Possible photospheric stellar contribution. 



with an age of ~ 120 Myrs, which is the median age of our sam- 
ple of galaxies, the predicted E'W values range from ~ 0.65 to 
~ 1.3 A (see iBrinchmann et al.ll2()08b !. Figure 2), depending on 
the metalUcity (0.4 < Z/Zq < 1.0). Since the metalUcity of the 
SF Gs in the GJVIASS sam ple has been estimated to be Z ~ 0.3Zq 
by iHallidav et al.l (|2008|) . our measured EW is fairly consistent 
with the model predictions. 

IVIost recent spectral population synthe sis codes include also 
the modeling of massive binaries system (lEldridge & Stanwavl 
^009; Brinchmann et al. 2008a), which may significantly effect 
the evolution of the WR stars. Up to now these models have only 
be tested in the local Universe, but their success in reproducing 



the observed spectral features also looks promising with regard 
to their apphcation to higher redshift samples. 

3.4.3. Fe II* 

In our composite spectrum, we were also able to identify some 
iron fine-structure emission lines, between 2000A and 3000A: 
Fe II* ^2365 A, Fe II* ^2396 A, Fe II* ^2612 A, and Fe II* 
/12626 A. These lines have been proposed to be generated by 
photon scattering in large-scale outflows (Rubin e t al.. .2010ai) . 
therefore their detection can be considered as a confirmation of 
the presence of outflows. 
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Fig. 6. Composite spectrum of the 74 spectra of the SFG sample (see Figure |5]l: zoom-in of the regions aroud the C IV and the Si 
IV doublets. The spectrum has been shifted into its stellar rest-frame. The vertical lines indicate the vacuum wavelength of the line. 
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Fig. 7. Average spectrum of the six galaxies with Lya line in emission. The spectrum has been shifted into the rest-frame given by 
the Lya. The expected position of the most prominent spectral lines is indicated. See Figure|5]for the colour legend. 



4. The UV continuum slope and dust extinction 

Under the assumption, which is quite robust when considering 
galaxies dominated by a young stellar population, that the shape 
of the UV continuum can be fairly accurately approximated by 
a power law 

Fa « Af^, 

where is the observed flux {erg i"' cm^^ A"') and /? is 
the continuum slope, Calzetti et al. (1994) showed that there is 
a strong correlation between the slope /? and nebular extinction 
measured in the optical using the Balmer decrement. 

For each of the 74 spectra in our sample of high-quality 
SFGs, the continuum was fitted with a least squares minimiza- 
tion, in the logiF^) - logiA) plane, of the average fluxes mea- 
sured in ten wavelength windows in the range between ~ 1250 A 
and 2600 A, as defined bv lCalzetti et al.l ([1994). These windows 
are chosen to avoid the strongest spectral features. In Figure |9] 



we present a spectrum from our sample, as an example. The 
range covered by the fitting windows is indicated by the vertical 
dotted lines in the upper panel, while the final fit is superimposed 
on the spectrum in red. The two lower panels show zoomed-in 
regions of the same spectrum, where the ten fitting windows (in 
the F,i - A plane) are indicated. The UV continuum slope (J3) 
distribution over the sample can be seen in Figure [10] The mean 
value for the sample is < yS >= -1.11 + 0.44 (r.m.s.). 

The definition of yS is not unique across the literature. First 
of all, it can be calculated dire ctly on spectra, as we did, or de- 
termined from the photometry (iBouwens et al.ll2009t) . Then, the 
range over which /3 is defined also changes from one work to 
another: when defined over a shorter wavelength range than that 
defined bv .Calz e tti et al. (1994), the quoted slope is systemati- 
cally redder rsee lCalzettfcoOlh . 

Keeping all this in mind, we compared the results of our 
analysis with others, for both local and high redshift galaxies, 
taken from the literature. The mean value of /3 over our sample is 
in good agreement with the < y6 > = -1.13 + 0.12 measuredby 
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Fig. 8. UV continuum slope as a function of the ratio between 
UV luminosity (not ex tinction c orrected) and IR luminosity. The 
red line shows the Meurer et al ] (11999,) relation. 



iHeckman et al.l (Il998l) from the spectra of 45 local star-forming 
galaxies observed with t he lU E satellite. For high redshift 
galaxies, iNoll & Pierinil (l2005h analyzed the spectra of 34 
star-forming galaxies at 2 < z < 2.5 from the FORS Deep Field 
(FDF) spectroscopic survey dNoU et al .1 12004*). The determined 
mean value of the continuum slope for the sub-sample of galax- 
ies with bluer spectra is < y6 > - -1.01 + 0.11, with which 
the mean value over our sample is fairly consistent, taking into 
account the different wavelength baseline (~ 1200 - 1800 A 
in their analysis, comp a red to ~ 1200 - 2600 A in this work). 
Finally, iBouwens et al.l 120091) presented a systematic measure- 
ment of the UV continuum slope over a wide range in redshift, 
using photometry. For galaxies at redshift z ~ 2, they find 
< P > ~ -1.5, with a scatter of ~ 0.30 (see iBouwens et ali 
120091 Figure 5), with which our results agree reasonably well, 
again considering the differences in methods and wavelength 
baselines used. 

iMeurer et al. l (fT999t) found an empirical relation between the 
ultraviolet spectral slope and dust extinction at 1600 A, in local 
star-forming galaxies given by 

Aisoo =4.43 + 1. 99 yS. 

We used this relation to estimate dust extinction for the 

galaxies of our sample . 

The lMeurer et al.1 (Il999l) relation was derived from an ob- 
served correlation between the ratio of infra-red luminosity over 
ultra-violet luminosity (not corrected for extinction) and the con- 
tinuum slope. We compared the continuum slopes computed 
for our sample with the ratio between luminosities measured 
at 1500A (see next section) and the IR luminosit ies extrapo- 
lated f rom the Spitzer-MlPS 24 /im fluxes (following'Pad di et al.l 
(|2^7!)), and o ur galaxies we r e foun d to be in reasonable agree- 
ment with the iMeurer et al. 1 119991) relation, as it is shown in 
Figure |8] 



Using the A leop values derived from the continuum slope, 
and a Calzetti law ( Calzetti et al ] 120001) . the color excess was 
derived. This determination of E(B - V) was then compared 
with other estimates: the extinction obtained by SED fitting 
an d the values c alcula ted from B and z bands using the relation 
of iDaddi et al.l (|2004|) . proposed to be valid for star-forming 
galaxies in the redshift range 1.4 < z < 2.5 

E(B-V) = 0.25 (B - z + 0.1), 

where B and z are AB magnitudes. 

Figure (TTj compares the different estimates. The color ex- 
cess from the continuum slope is quite consistent with the esti- 
mate derived from B and z bands, while it shows a larger disper- 
tion when compared to the estimate from SED fitting. However, 
in both cases the probability that the two compared data sets 
are linearly uncorrected, as derived by the Pearson coefficients 
(r^y = 0.57 and r^,. = 0.42), is < 0.1%. 

5. Star-formation rate 

After estimating the dust extinction from the continuum slope 
of the spectra in our galaxy sample, the star-formation rates 
from the rest-fra me UV luminosi ty were derived, by applying 
the relation from lKennicuti (119981) scaled for a Kroupa IMF 

SFRuv (Mo yr-i) = 0.9 x lO-^^ x L1500 A (^'■8 Hz '). 

First, the spectra were corrected for slit losses, so that in- 
tegrated fluxes calculated for specific bandpasses would 
match catalogue magnitudes obtained by photometry. The 
photometric values were found to be on average, about one 
magnitude brighter than the ones inferred from spectra. We 
checked whether a correlation existed between Amag (defined 
as magphnt - magspec) and the photometric magnitude, but 
found none. No correlation was found between Amag and the 
morphology of the objects. The absence of correlations between 
the magnitude corrections and galaxy properties likely indicates 
that the main cause of light loss in the spectral data is the slit 
width not encompassing the full seeing disk. 

For each of the corrected spectra, the 1500 A luminosity was 
obtained by measuring the average flux in a wavelength window 
ranging from 1350 A to 1650 A. Unextincted fluxes were then 
calculated using, for each galaxy, the estimate of dust extinction 
previously derived from its continuum slope. 

It is worth noting that wRy value had to be chosen to con- 
vert the esti mated extinction in to color excess, and we used the 
one from (ICalzetti et al.l 12000). But to correct the UV fluxes 
we directly used the Aieoo estimated from the UV continuum 
slope. Therefore, our extinction correction does not depend on 
the choice of a specific extinction law. 

The distribution of SFRs estimated from spectra, along with 
the ones computed by SED fitting, is presented in Figure [12] 
Throughout this paper, as pointed out in section IZTI all the SED 
fitting derived quantities were computed using Ma rastonI (l2005l) 
models with a Kroupa IMF. These models were chosen because 
they include a treatment of the thermally pulsing asymptotic gi- 
ant branch (TP-AGB) phase of stellar evol ution, but the SED 
fitting procedure was also performed using iBruzual & Chariot 
(120031) models with a IChabried (l2003l) IMF (see ICimatti et al. 
'2008"). We compared the SFR obtained from UV luminosity with 
both the estim ates from SED fitting, and we found the values 
computed with iMarastonI (|2005|) models to be the more consis- 
tent with the spectra-based estimate. 
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Fig. 9. Top panel: a spectrum from our sample with correspond- 
ing fit of the continuum superimposed in red. The vertical dashed 
lines delineate the range covered by the windows defined by 
ICalzetti et alj ( 1199 4) and used in the fitting procedure. Bottom 
panels: zoom over the wavelength range used in the fitting rou- 
tine, with the ten fitting-windows indicated. 



The galaxies of our sample show average < SFR > - 
52 Mq yr"', with r.m.s. - 48 Mp yr"', which is similar to the 
findings of Daddi et al] (l2007h for th eir samp le of K-sele cted 
star-forming galaxies at z ~ 2, and o f IShanlev et al.l (l20Q3h for 
their sample of LBGs at z ~ 3. 

5.1. The SFR-stellar mass correlation 

The existence of a correlation between stellar mass and SFR was 
also investigated, and the results are shown in Figure[T3] A posi- 
tive correlation between SFR and stellar mass has been reported 
to hold at different redshifts fe lbaz et al. 2007; Noeske et al. 
l2007HDaddi et al.l2007l:[s"antini et al..2009Lf annella et al..2009i: 




Fig. 10. UV continuum slope distribution over the sample of 74 
SFGs. The vertical red lines indicate the mean value for the sam- 
ple: </3>- -1.11, with r.m.s. - 0.44. 



iRodighiero et al.ll2010l) . despite the different techniques of SFR 
estimates used (radio luminosity and UV luminosity alone or 
combined with L/r), but the debate about its general validity 
and implications remains open. For our sample, we also found 
that more massive galaxies have higher SFRs (as estimated from 
LjjQQ ^), in agreement with previous studies. The big panel of 
Figure [13] shows the stellar mass-SFR correlation for our spec- 
troscopic sample of galaxies and for the photometric sample 



from D addi et al. ([2007) (for the sa ke of consistency, we scaled 



the values from lOaddi et al. I dloOTi) for our chosen IMF, before 
comparing the results). 

In order to estimate the slope of the SFR-Mass relation, 
a regression technique which treats symmetrically the vari- 
ables is required. In a case where the intrinsic scatter is likely 
to dominate any errors arising from the measurement pro- 
cess, the unweighted ordinary least-squares (OLS) bisector has 
been de monstrated to be a rec ommended choice (^Is obe et_d] 
1990; Fei gelson & Babu 1992; F eieelson & Jogesh Babull20nr 
In this method there is however the implicit assumption that the 
errors in the two quantities must be of the same order of magni- 
tude. In our study, we derived statistical errors on the SFR esti- 
mated from the spectra, and derived from simulations the errors 
on masses estimated using SED fitting techniques, and we found 
that both errors are on average about 30%, therefore the OLS 
bisector regression method is applicable. In Figure [T3]the black 
line, that comes from the linear fit to the data (on a logarithmic 
scale) of the GMASS sample, has a slope of 1.10 + 0.10, which 
is in agreem ent with the relation s found for other z ~ 2 samples, 
in particular lDaddi et al.l (12007'). who quote a slope of ~ 0.9 us- 
ing UV-derived SFRs, and |Pannella et al. (2009), who report a 
slope 0.95 + 0.07 and whose SFRs were computed from radio 
luminosity. We note that a similar proportionality is also seen at 
lower redshift (lElbaz et al.ll2007l: iNoeske et al.1l2007h . 

We also found that, in our sample, the SSFR decreases with 
the stellar mas s (see Figur e [13] insert panel), a trend with which 
the data from IDaddi et al.l (j2007) also seem to be consistent. A 
SSFR that declines with increasing stellar mass is found also in 
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Fig. 11. Comparison between E(B-V) derived from spectra UV 
continuum slope and, respectively, E(B-V) calculated from B 
and z magnitudes (top panel) and E(B-V) from SED fitting (bot- 
tom panel), with the one-to-one line in red. The Pearson correla- 
tion coefficient r^j, is also reported. 



iRodis rhiero et al ] (120101) for z ~ 2 sources, while IPannella et aP 
( 120091) and lDunne et all (l2009h report in their analyses that the 
SSFR of z ~ 2 star-forming galaxies is almost independent of 
stellar mass. 

However, we recall that several selection effects and biases 
might affect our sample of 74 SFGs. A more extended analysis 
of the SFR at z ~ 2 will be the presented in a forthcoming paper, 
for a larger sample of galaxies and using IR data provided by 
Herschel. 



Fig. 12. Distributions of extinction-corrected SFR [Mq yr '] de- 
rived from UV luminosity measured on spectra and extinction- 
ejected SFR [Mq '] from SED fitting (using .Marastonl 
(l2005h models with a Kroupa IMF). 



tral lines that originate from the absorption of the background 
stellar radiation by the interstellar medium. The most accurate 
estimate of the systemic rest-frame velocity of a galaxy is the 
one determined from spectral lines that are of stellar origin. The 
assumption that stellar spectral features are at rest with respect to 
their host galaxy is justified by considering that it is not possible 
to resolve stellar populations in high-redshift galaxies, and that 
the effects of stellar proper motions on a spectral line shift tend 
to statistically cancel out. The stellar lines observed in a galaxy 
UV rest-frame are extremely weak, and in the single GMASS 
spectra are almost undetectable because of the low S/N of the 
stellar continuum. It is possible to improve the S/N to identify 
the stellar spectral lines by working on composite spectra (see 
Figure [T4li. 

In Figure |5] the composite spectrum of 74 high-quality star- 
forming galaxies is shown, where two stellar and two nebu- 
lar lines can be detected and reliably measured: the absorption 
hnes C 111 ill 76 A, O IV ^1343 A, and the semi-forbidden C 
III] AA19Q6, 1908 A doublet and CII] ^2326 A emission Hne. 
Since the redshifts of the single galaxies are determined from 
the strongest lines, which are the interstellar absorption ones, the 
composite spectrum does not appear to be at rest with respect to 
the stellar and nebular features (see FigurefTSll. Therefore, in or- 
der to compute the velocity offsets of the lines tracers of the IS 
gas, first of all the composite spectrum was shifted into its rest- 
frame using the mean shift measured for the two chosen photo- 
spheric stellar lines and the two nebular carbon emission lines, 
weighted for their uncertainties (see Section l6^ . of ~ 135 + 22 
km/s.B 



6. Outflows 

We search for galactic-scale outflows in galaxy UV spectra by 
measuring the wavelength shift, of the centroid position, of spec- 



' The error in a weighted mean is defined as tr^ = 1/ where 
cr, are the measurement uncertainties associated with the stellar and 
nebular lines. 
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Fig. 13. Stellar mass as a function of SFR (big panel) for the GMASS sample analyzed in this paper (black dots). Stellar masses 
have been computed by means of SED fitting (using Maraston (2005) models with a Kroupa IMF). The extinction-corrected SFR 
[Mq yr^^] was derived from UV luminosity measured from spectra. The black line is a linear fit to the GMASS data (in logarithmic 
scale). We als o indicate the Pear son correlation coefficient r„.. The open red triangles represent the sample of z ~ 2 star-forming 
galaxies from lDaddi et al. |(|2003)(SFRsand masses scaled for a Kroupa IMF). In the little panel, stellar mass as a function of SSFR 
is shown. 



6.1. Velocity offset measurements 

The different kinds of lines observed in the UV rest-frame spec- 
trum of a SFG were already described in a previous section. 
Here we recall that, in the 1200 - 2900 A range, seven strong 
low-ionization interstellar absorption lines, tracers of the large- 
scale winds, are detectable, one of which is a blend; three low- 
ionization doublets are also detectable (Fe II and Mg II), in ad- 
dition to three high-ionization ones (Si IV, C IV, and Al III). 
The high-ionization doublets do not provide robust tools in the 
search for outflows, since at our spectral resolution an attempt to 
separate the ISM absorption from other components, such as the 
stellar wind P-Cygni profile, is impracticable. 

Therefore, in our analysis we concentrated only on the low- 
ionization interstellar absorption lines and doublets; in particu- 
lar: Si II /il260 A, C II ^1334 A, Si II ^1526 A, Fe II /11608 A, 
Al III ^1670 A, and Fe II i2344 A, and the doublets Fe II 
/12374, 2382 A, Fe II i2586, 2600 A, and Mg II ^^2796, 2803 A. 
The O H- Si II /II 303 A feature was excluded because it was im- 
possible to correctly determine the position of the centroids of 
the two blended lines. The other cited interstellar lines were all 
detectable and strong enough to ensure a quite good determina- 
tion of their centroid. 



Fe II /1 1608 A centroid is probably the least accurate, since 
it is the weakest of all the IS lines and its profile is broader 
and more asymmetric, these all being factors that may alter the 
position of the centroid during the fitting procedure. Moreover, 
we suspect that a blend of absorption features (Fe III, Al III, 
N II) typical of O and B stars in the 1600-1630 A range (see 
iKinnev et al.ll 1993b could play a major role in distorting the Fe 
II /11608 A line profile. 



6.2. Uncertainties 

For each line, two possible error sources related to the centroid 
measurement were estimated. The error due to spectral noise was 
computed by creating an artificial template with a continuum 
slope similar to that of the real composite spectrum and perfect 
gaussian spectral lines of similar intensities to those in the com- 
posite. Random noise was then added to the template until the 
S/N of the real composite spectrum was matched, and the cen- 
troids of the spectral lines of interest were measured. For each 
line, the r.m.s. of thirty repeated measurements, with respect to 
the input value, was then adopted as the measurement uncer- 
tainty. For the strongest absorption lines, the order of magnitude 
of this uncertainty is ~ 15-20 km/s; weaker lines have a higher 
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Fig. 15. Composite spectrum of the 74 spectra of the SFG sample (see Figure |5]i: zoomed image of the three detectable stellar 
photospheric absorption lines (C III /ill76 A, O IV /11343 A, S V /11501 A. The latter was not used in the analysis since it is too 
weak to be reliable measured), the C 111] AA1906, 1908 A nebular emission line doublet and the Cll] A2326 A nebular line. The 
vertical lines indicate the vacuum wavelength of the line. 



measurement uncertainty with a maximum of ~ 100 km/s for the 
O IV ^1343 A. 

We then evaluated how slight variations in the number of 
constituents of the composite spectrum would affect the aspect 
of the spectral lines or the position of their centroid, by creating 
about a dozen composite spectra and omitting from the proce- 
dure a number of single spectra randomly chosen between 5 and 
10. The lines of interest were then measured in all the spectra: 
the computed r.m.s. of the measurements for each line was even- 
tually defined as the "galaxy sampling uncertainty" and added 
quadratically to the aforementioned measurement error. On av- 
erage, the galaxy sampling uncertainty is ~ 10 km/s for the in- 
terstellar and nebular lines and ~ 35 km/s for the photospheric 
stellar ones. 

Finally, for unresolved C 111] ^^1906, 1908 A doublet, as al- 
ready mentioned, we fixed the intensities ratio of the two com- 
ponents to be 1.5, but we added to the line measurement an error 
component of ~ 39 km/s to take into account the un certainty in 
the r eal value of flux intensity ratio. For example, Ouider et al. 
(12009) measured a flux ratio F(1906)/F(1908) = 1.1 + 0.2 in 
the spectrum of the "Cosmic Horshoe", a gravitationally lensed 
star-forming galaxy at z ~ 2, whose physical properties such as 
SFR, mass, and dust extiction are similar to the average ones 
computed for our sample. We note that using a flux ratio of 1 . 1 
instead of 1.5, would not change our results. 

The total uncertainty for each line is quoted in Table |2] 



7. Results 

Figure[T6]shows zoomed-in spectral regions containing the main 
low-ionization interstellar absorption lines whose offsets were 
measured. The spectrum has been shifted into its rest-frame de- 
fined by the velocity correction of ~ 135 km/s computed from 
the stellar and nebular lines. 

Most of the lines are blueshifted with respect to their ex- 
pected positions, as indicated by the vertical dashed lines: this 
confirms the presence of outflows in the galaxies of our sample, 
at least on average. These results, qualitatively shown in the fig- 
ures, are presented in Table |2l A mean velocity offset of ~ 100 
km/s was computed for all the IS lines in our composite. The 
shift velocities measured for the single low-ionization lines in 
the composite spect rum are listed, along with the results from 
Sha plev et alJ (l2003h for the composite spectrum of a sample of 
Lyman-break selected galaxies (LBG) at z ~ 3. The EWs of all 
the lines were also measured, and presented. It was impossible 
to compare offsets and EWs of all the lines, since our composite 
spectrum extends more to the red than the LBG composite. 

In the wavelength range sampled by both spectra, given the 
uncertainties, the two sets of velocities are fairly consistent, 
probing that gas outflows are a common phenomenon in SFGs 
at all redshifts. 

On the other hand, the EWs measured in our composite are 
larger than in the LBG composite. Looking at the colour excess 
in the two samples, the LBG sample has E(B - V) values that 
are, on average, ~ 0.1 mag lower than the ones computed for 
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Fig. 16. Composite spectrum of the 74 spectra of the SFG sample (see Figure|5]): zoom-in of the strongest low-ionization interstellar 
absorption lines. The C IV high-ionization doublet is also labeled. The spectrum has been shifted into its rest-frame using the 
velocity correction derived from the stellar/nebular lines (~ 135 km/s). The vertical lines indicate the expected positions of the line. 



our sample. This can be intepreted as evidence of the correlation 
between the interstellar absorption line strength and dust extinc- 
tion. 



7.1. Searching for correlations between UV spectroscopic 
properties and galaxy physical properties 

The spectral analysis of the composite spectrum has provided 
average information about our entire sample of z ~ 2 high- 
quality star-forming galaxies. The next step was to investigate 
possible correlations between galaxy UV spectral characteris- 
tics (EW, outflows) and galaxy general properties, in particular 
stellar mass, SFR, and dust-extinction. 

Stellar masses were computed by means of the SED fitting 
of iMarastonI ( 120051) models, while the colour excess was derived 



from the continuum slope and the SFR was obtained from UV 
spectroscopic luminosities, as described in Section|5] 

Our sample was divided into two equally populated bins, for 
each property, and a composite spectrum of each group of galax- 
ies was created, following the same procedure adopted in the cre- 
ation of the composite of the total sample. We decided to work 
with only two bins because of the relatively small number of 
spectra at our disposal. As already stated, to have reliable mea- 
surements of centroids and EWs of the spectral lines of interest, 
a high S/N is needed in the composite spectra, and the compo- 
sition procedure is the more effective in improving the S/N the 
higher the number of single spectra involved. 

Figure[T7]shows the distribution of the chosen physical prop- 
erties over our sample, with the binning division values marked 
and the median values of each bin reported. The composite spec- 
trum of each bin was shifted into the rest-frame given by the 
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Table 2. Velocity offsets and EWs of strong low-ionization lines, measured in the composite spectrum of the full SFG sample. 



Ion 


Vacuum 


IE (eV)" 


EW (A)'" 


V (km/s) 




wavelength (A) 




GMASS LBG' 


GMASS 


LBG' 


Sill 


1260.42 


8.15 


1.47 ±0.09 1.63 ±0.10 


-79 ± 19 


-110 ± 60" 


O I - Sill'' 


1302.17, 1304.37 


- 


2.43 ±0.26 2.20 ±0.12 


-123 ±42 


-270 ± 60'' 


CII 


1334.53 


11.26 


2.09 ±0.15 1.72 ±0.11 


-105 ± 33 


-150 ± 60' 


Sill 


1526.71 


8.15 


1.94 ±0.08 1.72 ±0.18 


-79 ± 19 


-110 ±60' 


Fe II 


1608.45 


7.87 


1.31 ±0.11 0.91 ±0.15 


24 ±40 


-60 ± 60" 


Al II 


1670.79 


5.99 


1.54 ±0.08 1.04 ±0.15 


-78 ± 27 


-100 ±60' 


Fell 


2344.21 


7.87 


2.45 ±0.13 


-139 ±27 




Fell^ 


2374.46 


7.87 


1.87 ±0.13 


-72 ± 32 




Fell^ 


2382.76 


7.87 


2.21 ± 0.15 


-72 ± 27 




FelP 


2586.65 


7.87 


2.65 ± 0.17 


-108 ± 27 




FelF 


2600.17 


7.87 


2.77 ± 0.14 


-108 ± 28 




Mg II'' 


2796.35 


7.65 


2.43 ± 0.10 


-176 ± 27 




Mg II'' 


2803.53 


7.65 


2.07 ± 0.10 


-176 ± 27 





Notes. Ionization energy. 

Rest-frame equivalent wid t hs. 

Values from Shap lev et all j2003h . 
*''' All the mea su remen ts, both ours and Shapley's, refer to the blend of the two lines. 
"^'' IShaplev et alj (l2003h quotes only the mean error in the velocity offsets. 
(/) te) (*) Doublets. 




Fig. 17. Distribution over the sample of stellar mass, SFR (from 1500 A luminosity measured on the spectra with dust extinction 
derived from the continuum slope) and E(B-V) (computed from the UV continuum slope). The median value of both the sub-samples 
is also reported in each panel. 



stellar and nebular lines. Then, the velocity offsets of the low- 
ionization absorption lines were measured, which are presented 
in Table[3] Since Table[3]is a bit "dense", the mean velocity off- 
set and EW of the interstellar low-ionization absorption lines in 
each bin are reported in Tabled 

7.1.1. Outflow velocities 

The IS lines show a significant velocity shift in all the composite 
spectra {6vis ~ 100 km/s), but there seem to be no evident corre- 
lations between the velocity offsets of the interstellar absorption 
lines and the physical properties of the galaxies in our sample, 
apart from some weak differences in the colour excess and SFR. 

However, our results do not completely exclude the exis- 
tence of these correlations, since we are aware of some factors 
that may hide them. First of all, the magnitude of the uncer- 
tainties might be playing a significant role in hiding real weak 
correlations. Moreover, the distributions of SFR, stellar mass. 



and colour excess do not favor a division of the sample into 
two equally populated bins, since most of the galaxies are con- 
centrated close to the threshold value, and the eiTors are large 
enough to cause confusion between the two bins. Finally, the 
fact has to be considered that there may be absorption also from 
gas which is not out flowing, expecially in more massive galaxies 
(ISteidel et al.ll20ldl) . This v ~ component may affect the emer- 
gent line profile in the sense that the velocities computed from 
the centroid shift of the IS lines may be underestimated. 



7.1 .2. Equivalent widths of outflow-related absorption lines 

The EW measurements show significant differences, when com- 
paring the composite spectra of different sub-samples. A positive 
trend between the EWs of the interstellar absorption lines and 
SFR, stellar mass, and E(B-V) is suggested, which is the main 
result of this last part of our analysis. 
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Table 3. Equivalent widths and velocity offsets of the strongest interstellar absorption Unes (plus the He II /11640 A emission line), 
measured in composite spectra. 



Ion 


Vacuum 


EW 


V 


EW 


V (km/s) 




A h 

Av 




wavelength (A) 


(A) 


(km/s) 


(A) 


(km/s) 


(A) 


(km/s) 




SFR 


BinP: 20.24 M^yr-' 


Bin2'^: 61.06 M^yr"' 






Sill 


1260.42 


A m 1 A AA 


m 1 1 A 


1 C7 1 A AA 

1.5/ ± u.uy 


1 /I C 1 1 A 

-14j ± 19 


1.62 ± 0.33 


-49 ± 27 


CII 


1334.53 


1.98 ±0.15 


-114 ±33 


2.17 ±0.15 


-138 ± 33 


1.10 ±0.33 


-24 ± 47 


Sill 


1526.71 


2.44 ± 0.08 


-97 ± 19 


3.73 ± 0.08 


-145 ± 19 


1.53 ± 0.27 


-49 ± 27 


Fell 


1608.45 


1.09 ± 0.11 


9 ± 40 


1.21 ±0.11 


-47 ± 40 


1.11 ± 0.25 


-56 ± 57 


Hell 


1640.42 


-1.13 ± 0.10 




-0.87 ± 0.10 




0.77 ± 0.16 




AlII 


1670.79 


1.35 ± 0.08 


-108 ± 27 


1.77 ±0.08 


-95 ± 27 


1.31 ± 0.22 


13 ±38 


Fell 


2344.21 


2.19 ± 0.13 


-185 ± 27 


3.02 ±0.13 


-120 ± 27 


1.38 ±0.38 


65 ±38 


Fen 


2374.46 


1.81+0.13 


-150 ± 32 


2.77 ±0.13 


-89 ± 32 


1.53 ± 0.43 


61 ±45 


Fen 


2382.76 


2.06 ± 0.15 


-150 ± 27 


2.28 ± 0.15 


-89 ± 27 


1.11 ±0.33 


61 ±38 


Fen 


2586.65 


1.88 ±0.17 


-159 ± 27 


3.74 ± 0.17 


-87 ± 27 


1.99 ± 0.84 


72 ±38 


Fen 


2600.17 


2.57 ±0.14 


-159 ± 28 


3.33 ± 0.14 


-87 ± 28 


1.30 ±0.38 


72 ±40 


MgH 


2796.35 


2.53 ± 0.10 


-239 ± 27 


2.18 ±0.10 


-137 ± 27 


0.86 ± 0.17 


101 ± 38 


MgH 


2803.53 


1.87 ± 0.10 


-238 ± 27 


2.37 ±0.10 


-137 ± 27 


1.27 ± 0.26 


101 ± 38 


Stellar Mass 


Binr: 9.67 log(Mo) 


Bin 2": 10.04 log(Mo) 






Si II 


1260.42 


1 61 + 09 


-140 ± 19 


1 89 + 09 


-138 ± 19 


1.17 ±0.21 


2 ±27 


cn 


1334.53 


2.25 ± 0.15 


-134 ± 33 


234 + 0.15 


-91 ± 33 


1.04 ±0.31 


43 ±47 


Sin 


1526.71 


2.47 ± 0.08 


-140 ± 19 


3.14 ± 0.08 


-138 ± 19 


1.27 ± 0.21 


2 ±27 


Fen 


1608.45 


1.17 ±0.11 


46 ±40 


1.28 ±0.11 


-13 ± 40 


1.09 ± 0.24 


-59 ± 57 


Hen 


1640.42 


-1.23 ±0.10 




-0.81 ± 0.10 




0.66 ± 0.14 




Ain 


1670.79 


1.33 ± 0.08 


-80 ± 27 


1.87 ± 0.08 


-92 ± 27 


1.41 ± 0.24 


-12 ±38 


Fen 


2344.21 


2.37 ± 0.13 


-149 ± 27 


2.72 ±0.13 


-139 ±27 


1.15 ± 0.30 


9 ±38 


FeH 


2374.46 


1.68 ± 0.13 


-82 ± 32 


2.51 ±0.13 


-93 ± 32 


1.49 ± 0.42 


-11 ± 45 


FcH 


2382.76 


1.71 ± 0.15 


-82 ± 27 


2.90 ±0.15 


-93 ± 27 


1.7 ±0.58 


-11 ± 38 


FeH 


2586.65 


2.59 ± 0.17 


-139 ± 27 


3.31 ±0.17 


-101 ± 27 


1.28 ± 0.45 


37 ±38 


Fen 


2600.17 


2.88 ± 0.14 


-139 ± 28 


3.27 ±0.14 


-101 ± 28 


1.14 ±0.32 


37 ±40 


Mgn 


2796.35 


2.24 ±0.10 


-183 ± 27 


2.37 ± 0.10 


-209 ± 27 


1.06 ± 0.21 


-26 ± 38 


Mgn 


2803.53 


1.63 ±0.10 


-183 ± 27 


2.66 ± 0.10 


-209 ± 27 


1.63 ± 0.37 


-26 ± 38 




E(B-V) 


Bin Y : 


0.16 


Bin 2'^ : 


0.27 






Sin 


1260.42 


1.43 ± 0.09 


-84 ± 19 


2.25 ± 0.09 


-38 ± 19 


1.57 ± 0.31 


46 ± 27 


cn 


1334.53 


1.85 ± 0.15 


-160 ±33 


2.18 ±0.15 


-58 ± 33 


1.18 ± 0.36 


102 ± 47 


Sin 


1526.71 


1.70 ± 0.08 


-84 ± 19 


2.27 ± 0.08 


-38 ± 19 


1.34 ± 0.22 


46 ±27 


Fen 


1608.45 


1.25 ± 0.11 


7 ±40 


1.61 ±0.11 


86 ±40 


1.29 ± 0.29 


78 ±57 


Hen 


1640.42 


-0.88 ±0.10 




-0.96 ±0.10 




1.10 ± 0.22 




Ain 


1670.79 


1.31 ± 0.08 


-46 ± 27 


1.81 ±0.08 


-51 ± 27 


1.38 ± 0.23 


-5 ±38 


Fen 


2344.21 


2.28 ±0.13 


-103 ±27 


2.76 ± 0.13 


-110 ±27 


1.21 ± 0.32 


-7 ±38 


Fen 


2374.46 


1.59 ±0.13 


-78 ± 32 


2.15 ±0.13 


-21 ± 32 


1.35 ± 0.37 


57 ±45 


Fen 


2382.76 


1.96 ± 0.15 


-78 ± 27 


2.40 ±0.15 


-21 ± 27 


1.23 ±0.38 


57 ±38 


Fen 


2586.65 


1.53 ± 0.17 


-98 ± 27 


3.89 ±0.17 


-58 ± 27 


2.55 ± 1.28 


39 ±38 


Fen 


2600.17 


2.11 ± 0.14 


-98 ± 28 


3.58 ±0.14 


-58 ± 28 


1.70 ± 0.54 


39 ±40 


Mgn 


2796.35 


1.58 ±0.10 


-186 ± 27 


3.04 ±0.10 


-124 ± 27 


1.93 ± 0.47 


63 ±38 


Mgn 


2803.53 


1.64 ±0.10 


-186 ± 27 


2.32 ± 0.10 


-123 ± 27 


1.41 ± 0.30 


63 ±38 



Notes. EW ratio = EWmghBin I EWlowBin 
Av = VHighBin - VlowBin 

Median of each bin. 



Table 4. Mean measured equivalent widths and velocity offsets. 







EW„ 


(A) 






V75(km/s) 




FW * 

^ ratio 


Av' 




SFR(Moyr-') 


Bin 1" : 


20.24 


Bin 2" : 


61.06 


Bin 1" : 


20.24 


Bin 2" : 


61.06 










1.90 ± 


0.15 


2.51 ± 


0.24 


-141 : 


1 19 


-110 


±9 


1.33 ± 


0.16 


31 ± 


21 


Stellar Mass (log(Mo)) 


Bin l" : 


: 9.67 


Bin 2" : 


10.04 


Bin 1° ; 


: 9.67 


Bin 2" : 


10.04 












1.99 ± 


0.16 


2.52 ± 


0.18 


-117 : 


t 18 


-118 


± 16 


1.26 ± 


0.13 


-1 ± 


24 


E(B-V) 


Bin r : 


: 0.16 


Bin 2" : 


: 0.27 


Bin l" : 


: 0.16 


Bin 2" 


: 0.27 












1.69 ± 


0.09 


2.52 ± 


0.20 


-99 ± 16 


-51 ± 16 


1.50 ± 


0.14 


48 ± 


23 



Notes. Median of each bin. 

(c) A _ 

" ' - VHighBin - VLowBin 
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Fig, 14. Comparison of the total average spectrum (black, below; 
the same seen in Figure |5]l of our sample with the spectrum of 
a single galaxy (red, above). See Figure |5] for the colour code 
adopted. 

In other samples of star-forming galaxies, at different red- 
shifts but with comparable properties of stellar mass and SFR, 
interstellar absorption lines associated with large scale outflows 
have been found to be stronger in galaxies wit h higher SFR, 
highe r stellar mass, and more heavil y extin cted (IShaplev et alj 
l2003HWeiner et al .1120091: iRubin et aLlfe oiObV Given the correla- 
tions between these galaxy properties, it is not possible to estab- 
lish the characteristic with which inter-stellar lines parameters 
are more strongly connected, even if the most significant of the 
three correlations appears to be the one with E(B-V). 

From these results it may be inferred that the EWs of the 
absorption inter-stellar lines are probably related to the velocity 
dispersion of the gas, or to its geometry, or to a combination of 
the two effects, all of which are difficult to distinguish with the 
actual data. 

Finally, we report the EW of He 11 /II 640 A to be stronger 
in the composite spectra of the less massive and star-forming 
galaxies. 

8. Summary and conclusions 

We have studied a sample of 74 star-forming galaxies at z ~ 2 
from GMASS by analyzing their rest-frame UV spectra. The first 
part of this work focused on the determination of dust extinc- 
tion and SFR for the galaxies in our sample, while in the sec- 
ond part large-scale gas-outflow-related properties were also an- 
alyzed. The main results of these analyses may be summarized 
as follows. 

1 . The continuum slope of the 74 individual spe c tra of our sam- 
ple was measured, following ICalzetti et al.l d 19941) . with a 
mean value for the sample of < /3 >= -1.11 ± 0.44 (rm.s.); 
this value confirms the results of similar studies in the same 
redshift range of our sample, which are also consistent with 
the findings in the local Universe. 

2. Dust extinction was derived for each galaxy from the contin- 
uum slope of its spectrum, and used to correct the flux mea- 



sured at 1500 A (rest-frame). The corrected flu x was then 
converted into SFR by applying the relation from lKennicutll 
(!l998) scaled for a Kroupa IMF; the galaxies of our sample 
show average < SFR > = 52 + 48 Mq yr"' (rm.s.). 

3. A positive correlation between SFR and stellar mass was 
found, with a logarithmic slope of the relation of 1 . 10 + 0. 10, 
in agreement with previous studies; the SSFR was also found 
to decrease with stellar mass. 

4. A composite spectrum was created by averaging the high- 
quality spectra of 74 SFGs, to obtain average information 
about large-scale gas-outflow-related properties. The low- 
ionization absorption lines associated with the interstellar 
medium were found to be blueshifted, with respect to the rest 
frame of the system, indicating the presence of outflowing 
gas with velocities of the order of ~ 100 km/s. Our measured 
velocity offsets were found to be fairly consistent with that 
obtained by Shaple v et alj (i2003l) for the composite spectrum 
of a sample of LBGs at z ~ 3, in the wavelength range sam- 
pled by both spectra (between 1000 and 2000 A). 

5. The EWs of the low-ionization absorption lines associated 
with the interstellar medium were also measured: the EWs 
measured in our composite appear to be larger than in the 
LBG composite. This could be interpreted as evidence of 
a trend between the interstellar absorption-line strength and 
dust extinction, since the LBGs have E{B - V) values that 
are, on average, ~ 0. 1 mag lower than the ones computed for 
our sample. 

6. In the composite spectrum of the galaxies showing Lya 
in emission, a velocity difference of ~ 566 + 27km/s be- 
tween this line and the IS absorption ones was measured. 
Moreover, the comparison of composite spectra of two sub- 
samples, one with Lya in emission, the other with the fea- 
ture in absorption, showed that the continuum slope is bluer 
in the first one (respectively, PemLya — -1.40 ± 0.06 and 
PabsLya - -1.04 + 0.03), thus Confirming the existence of 
a relation between UV continuum slope, and therefore dust 
extinction, and the Lya profile in SFGs. 

7. Possible correlations between galaxy UV spectral character- 
istics (EW, outflows) and galaxy general properties, in partic- 
ular stellar mass, SFR, and dust extinction, were investigated 
by dividing the sample into bins, for each of the aforemen- 
tioned properties. No evidence of significant correlation was 
found between the velocity offsets of the interstellar absorp- 
tion lines and the physical properties of the galaxies in our 
sample. However, this does not completely exclude the exis- 
tence of such correlations, since the magnitude of the uncer- 
tainties may hide them. 

8. The He 11 /11640 A emission line, which is a feature of winds 
from WR stars, was measured in all the composite spectra 
and the values were compared with the predictions of the 
most recent stellar population synthesis models. Given the 
physical properties of the galaxies of our sample, the models 
predict EWs which are fairly consistent with our measure- 
ments. We also report the He 11 /11640 A to be stronger in 
the composite spectra of the less massive and star-forming 
galaxies, a puzzling result since the higher fraction of WR 
stars is expected at higher metallicities. 

9. A positive trend between the EWs of the interstellar absorp- 
tion lines and SFR, stellar mass, and E(B-V) has been found, 
which has been reported to also hold at different redshifts. 
The most significant of the three correlations appears to be 
the one with E(B-V), which suggests that the EWs of the ab- 
sorption inter-stellar lines are probably related to the velocity 
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dispersion of the gas, to its geometry, or to a combination of 

the two effects. 

New data provided by Herschel will allow deeper investiga- 
tions of the SFR in the high-redshift galaxies of the GMASS 
sample and, consequently, help us to constrain the evolution of 
the star formation over cosmic time. 

Moreover, to improve our understanding of the phenomenon 
of large-scale outflows, deeper spectroscopic surveys will be 
needed in the future. In particular, it would be extremely use- 
ful to have multi-wavelength spectral data. High-quality spec- 
tra sampUng not only the rest-frame UV range, but also longer 
wavelengths, could allow a robust determination of systemic ve- 
locities and, consequently, the study of the warm phase of large- 
scale outflows in single galaxies, a task that to date has only been 
attempted for lens-magnified objects. 
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